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Abstract. Electromigrated nanoscale junctions have proven very useful for 
studying electronic transport at the single-molecule scale. However, confirming 
that conduction is through precisely the molecule of interest and not some 
contaminant or metal nanoparticlc has remained a persistent challenge, typically 
requiring a statistical analysis of many devices. We review how transport 
mechanisms in both purely electronic and optical measurements can be used 
to infer information about the nanoscale junction configuration. The electronic 
response to optical excitation is particularly revealing. We briefly discuss surface- 
enhanced Raman spectroscopy on such junctions, and present new results showing 
that currents due to optical rectification can provide a means of estimating the 
local electric field at the junction due to illumination. 
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1. Introduction 

Electronic devices with molecules and organic components as active elements offer 
new limits of device scaling and functionality, and are also of fundamental physical 
interest. Studies at the single molecule level probe the physics of electronic conduction 
and optical interactions in regimes that have been previously inaccessible. Three- 
terminal electronic measurements have proven invaluable in many systems, enabling 
electronic transport to function as a spectroscopy of available states. [1] While much 
progress has been made using two-terminal measurements to examine charge transport 
through molecules [2-20] , three-terminal measurements greatly increase the available 
information in single-molecule measurements. For example, through the use of a gate 
electrode it is possible to study conduction as a function of the molecule's redox 
potential, analogous to what can be done in electrochemical experiments. 

A common approach to fabricating three-terminal single-molecule devices 
is known as the "electromigration technique" , in which thermally assisted 
electromigration is used to create a nanoscale gap in a narrow metallic wire situated 
on a gate insulator material [21]. If a molecule of interest resides in the nanogap, the 
broken ends of the wires are used as the source and drain contact electrodes, and an 
underlying conductive substrate can be used as a gate electrode. The resulting device 
is a single- molecule transistor (SMT). Electromigrated breakjunctions have been used 
to study transport through individual nanocrystals [21], small metal particles [22], a 
variety of molecules [23-39] , and even individual atoms [40] . 

A central hurdle in most single-molecule electronic measurements is to 
demonstrate unambiguously that transport is occurring through only a single molecule 
of interest. This is complicated by the lack of direct imaging techniques with sufficient 
resolution, except for the two-terminal example of scanning tunneling microscopy 
(STM), which allows imaging with atomic resolution on conducting substrates. In this 
paper we review SMT fabrication (Section[2]), and then discuss device characterization 
based on electronic transport alone as a diagnostic (Section [3]). In that circumstance 
it is necessary to search for features characteristic of conduction through a single 
molecule in addition to attributes uniquely identifying the molecule of interest. 
Generally the search for single molecule devices requires a statistical approach and 
many control experiments. Exciting recent work (Section [4J demonstrates that 
simultaneous single-molecule optical spectroscopy and transport is possible. In such 
measurements surface-enhanced Raman spectroscopy (SERS) can give the vibrational 
spectrum of the specific molecule through which transport is taking place. We present 
new data using nonlinear transport and optical rectification to arrive at a quantitative 
estimate for the local optical field experienced by the molecule under illumination. 

2. Device Fabrication 

Fabrication begins with the preparation of arrays of devices on an n + Si substrate 
with a 200 nm SiC>2 insulating layer. 1 nm of Ti and 15 nm of Au are evaporated onto 
nanowire patterns defined by electron beam lithography. Bowtie-shaped constriction 
patterns are produced with minimum widths of approximately 100 nm (Figure [T^,). 
Additional gold pads for contacting the source and drain electrodes arc also defined by 
electron beam lithography. The array of samples is cleaned in solvents then exposed to 
an oxygen plasma for 1 min to remove trace organics. The molecules are dispersed for 
approximately monolayer coverage, or allowed to self-assemble onto the Au surface, 
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depending on the molecule in question. 

Elcctromigration (with some thermal enhancement due to Joule heating) is used 
to break the constriction into distinct source and drain electrodes, with the intent 
that the molecule will reside in the resulting nanoscale gap. This process has been 
studied extensively by a number of investigators [41-46] . At large current densities in 
the constriction, the electrons transfer sufficient momentum to the lattice to move 
atoms at the surface and at grain boundaries. As a constriction is formed the 
local current density increases, leading to a runaway migration and the formation 
of an interelectrode gap. Elevated local temperatures enhance the rate of this 
process through the increased diffusion of metal atoms. To mitigate concerns about 
adsorbed contaminants we prefer to electromigrate junctions at low temperatures in a 
cryopumped ultrahigh vacuum (UHV) environment, though this has not been possible 
in the optical experiments discussed in Section |4] The critical gap size, defined as 
the minimum separation between the two resulting electrodes, is typically 1-3 nm. 
With good, reproducible lithography to produce the constrictions and an appropriate 
elcctromigration procedure, yields of such small gaps can exceed 90%. 




Figure 1. (a) SEM image of a gold nanowire structure on Si/Si02 prior to 
electromigration. Inset shows a magnified SEM image of a similar nanowire after 
the breaking procedure, (b) Ideal device configuration after the electromigration 
process is complete for a Cgo-based junction. 

An assortment of voltage ramping techniques have been utilized by different 
research groups in an effort to form the smallest, cleanest, and most consistently- 
sized nanogaps [41-46]. All techniques involve minimizing the scries resistance to 
avoid overheating and subsequent melting of a nanowire; and an immediate reduction 
in applied bias across the break after the gap has opened in order to avoid creating 
an excessively large interelectrode spacing. We have found success with a feedback 
controlled method of repeatedly ramping the source-drain voltage, Vsd- 
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The desired final device configuration is shown schematically in Figurc[T|3. Subtle 
differences in junction configuration down to the Angstrom level can strongly influence 
device properties. The WKB approximation dictates that the current through a metal- 
vacuum-metal gap decreases by roughly a factor of e 2 (~7.4) for every Angstrom 
increase in the size of the vacuum gap [49] . This exponential dependence of tunneling 
on interelectrode distance ensures, as in an STM. that the tunneling electrons probe 
a volume containing at most one or two molecules. Precisely how molecules that 
begin on the surface of the constriction end up in the interelectrode gap is not 
known. The inherent randomness in the device formation process leads to many 
electrode-molecule configurations and resulting electrical conduction characteristics 
for any group of devices. Historically this has mandated a statistical approach to 
device characterization based on transport. By analysis of the transport data, as 
explained below, it is often possible to infer the number of molecules in the tunneling 
region and the relative couplings of the molecule to the source and drain electrodes. 
Details about molecular positioning and bonding remains unavailable, though new 
measurements such as those in Sect. [4] contain much information. 

3. Low temperature electronic transport 

The differential conductance, dl/dV, as a function of Vsd, calculated or measured 
using quasi-dc techniques, is the main transport tool to assess the nature of a 
given device. There are four main classes of differential conductance traces, with 
characteristic abundances determined by examining large numbers (thousands) of 
junctions [37] (though these vary depending on molecule type), f) Devices with no 
measurable differential conductance indicate an interelectrode breakjunction gap that 
is too large for measurable conduction by tunneling (~ 10%). 2) Devices with a weakly 
nonlinear conductance plot wherein dl/dV drops around Vsd = 0, but does not form 
a clear blockaded region denote brcakjunctions with small to moderately sized gaps, 
but without an active element within the nanogap (Figure [2^ ~ 45%). 3) Devices 
that show a region of zero conductance bordered by sharp peaks most likely have a 
gap containing a molecule or small particle (Figure (2b, percentage depends strongly 
on molecule type) . 4) Devices with a zero- bias resonance suggest a stronger molecule- 
electrode coupling (percentage depends strongly on molecule type, ranging from zero 
for alkane chains to several percent for transition metal complexes containing unpaired 
d electrons [31]). 

Further investigation is required to identify those devices that contain a single 
molecule of interest. Confounding possibilities include surface contaminants, metal 
grains produced during the electromigration process, or simply more than one molecule 
of interest. As explained below, the observed gate dependence of conduction and an 
understanding of the relevant transport physics has enabled progress in eliminating 
spurious devices. The final yield of gateable likely SMTs is typically 10-15% of the 
starting devices, based on a total sample size of thousands of devices. 

3.1. Coulomb blockade 

In analogy to conventional single-electron transistors, the SMT can be thought of 
as consisting of 6 parts: source, drain and gate electrodes, the molecule, and the 
tunneling connections to the source and drain. The tunneling barriers are established 
by the geometry and chemistry of the molecule/electrode interface. We assume 
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that the polarizability of the complex is much larger than that of the insulating 
barriers; therefore the voltage drops on the barriers are proportional to their respective 
coupling strengths Tg and . Elementary transport characteristics can be calculated 
by treating tunneling barriers as a capacitive and a resistive element [50]. The 
positive and negative slopes of the blockaded region (the tunneling thresholds) are 
determined by the ratios Cg/Cd and Cq/Cs, respectively [30] (Cg, Cd, and Cs 
are the capacitances between the molecule-gate, molecule-drain, and molecule-source 
electrodes). 

Traces of dl/dV measured as a function of Vsd, obtained within a range of 
incrementally changing gate voltage, Vq, can be compiled to form a map known as a 
stability diagram, illustrating where tunneling is both allowed and prohibited. Not all 
devices exhibiting a zero-bias resonance or a conductance gap will show dependence 
on applied gate bias, Vq- This poor gate coupling may occur as a result of the detailed 
device geometry or possibly as a result of the orientation of the active element within 
the brcakjunction gap [31,51]. Stability diagrams of devices that do exhibit gate 
dependence are a primary means of displaying pertinent transport data reflecting 
important aspects of SMT behavior. 

In equilibrium the source and drain electrodes have a common electrochemical 
potential, /x, the Fermi level. When the source and drain electrodes are biased such 
that the source electrode has a chemical potential /ii and the drain electrode has 
chemical potential /i 2 , current will flow as long as the Coulomb charging energy (E c ) 
and the discrete energy level spacing (A) have been overcome, and an energy level 
of the molecule lies between [i\ and ^2 (at finite temp, this range is extended by the 
thermal energy ifceT [52]). This current flow is due to resonant tunneling of electrons 
from the source to the lowest available single particle state producing a dl/dV peak 
at the edge of the conductance gap (Figure [3];). Gate voltage shifts the active element 
levels with respect to [i\ and /i2- The resulting stability diagram is determined by 
both the active element's spectrum and its capacitive couplings to the source, drain, 
and gate electrodes. 

When the active element has a negligible single-particle level spacing (e.g., a few- 
nm metal grain), the discrete spectrum seen in transport arises purely from E c , the 
Coulomb repulsion that must be overcome to add another electron to the element. 
The result is a series of regularly spaced diamonds in the stability diagram, and the 
suppressed conduction in the absence of a discrete level available for transport is 
called Coulomb blockade. Each diamond region represents a different charge state 
of the particle (i.e. a stable region of fixed average electron number). In this case 
the electron addition energy depends upon the classical capacitance of the particle, 
including interactions with the electrodes. 

Frequently the single-particle electronic level spacing, A, cannot be neglected 
when considering transport. If the island is a molecule such as Ceo, the discrete energy 
level spacing between single particle energy states at equilibrium occurs between the 
highest occupied molecular orbital level (HOMO level) and the lowest unoccupied 
molecule orbital level (LUMO level), commonly called the HOMO-LUMO gap. We 
can see from equation (1) that the size of the Coulomb diamond centered around 
^sd, Vg = depends on both the Coulomb charging energy as well as the magnitude 
of the HOMO-LUMO gap. 




(1) 
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where C cq = Cq + C$ + Co, dVsn/dVG is the slope of the diamonds, and AVq is the 
voltage spacing between successive charge states. 

where A is the spacing between two discrete energy levels (e.g. the HOMO-LUMO 
gap), and e 2 /C eq is the bare Coulomb charging energy, E c . 

Unlike metal particles, most molecules have a very limited number of accessible 
charge states, which do not typically appear with even spacing. Solution-based 
electrochemistry (e.g. cyclic voltammetry) demonstrates that molecules have a limited 
number of accessible valence states. [37] If one finds that it is possible to add many 
electrons to a potential device via gating, or if the spacing between charge state 
transitions is highly regular, it is extremely unlikely that the active region of the 
potential device contains a single small molecule. 

Stability diagrams with linear tunneling barrier edges possessing only two 
characteristic slopes (one positive and one negative) can be simulated using the basic 
capacitance model of sequential electron tunneling [36] . A system with more that one 
active element will always have more than two characteristic slopes. This provides a 
means of assessing whether tunneling is occurring through multiple molecules. 

When scrutinizing a potential SMT with a clear blockaded region, one must also 
look at the electron addition energy (the max source-drain bias of the blockaded 
region), which should be of a sensible size for the molecule in question. The classical 
capacitance of a molecule is so small that an electrostatic charging energy over 100 meV 
is not unreasonable, even without taking into account molecular level spacing. 

Some characteristics in stability diagrams may uniquely identify an SMT made 
with a particular molecule. At sufficiently high source drain bias, an electron may 
tunnel from the source into an excited single particle state of the island (Figure [3)i) . 
Additional dl/dV peaks that parallel the edges of the blockaded diamonds in stability 
diagrams appear beyond the edge of the conductance gap when tunneling occurs 
through newly accessible quantized excitations. The location of each dl/dV peak 
outside the conductance gap provides information on the excitation spectrum of a 
SMT. 

The conductance gap disappears (ordinarily) at the charge degeneracy point, 
Vg = V c , where the total energy of the system is the same for two different charge 
states of the molecule. When the gate voltage passes V c as it is ramped in the positive 
direction, the average number of charges on the molecule changes by one electron. The 
absolute equilibrium electron population, n, is determined by Vq, molecule/electrode 
intcrfacial charge transfer, and the local charge distribution around the molecule. 
Unfortunately, n cannot be determined solely from the data obtained in the stability 
diagram. 

Each dl/dV peak on the Vq < V c side of a stability diagram represents an opening 
of a new tunneling pathway where an electron tunnels onto the n-electron state, for 
example, to transiently generate the (n + l)-electron ground or excited states. In this 
way, the peaks probe the excitation energies of the molecule. Conversely, each dl/dV 
peak that appears at Vq > V c probes the ground and excited states of the n-electron 
molecule [24] . The energy of these quantized excitations can be determined from the 
source/drain voltage at which they intercept the conductance gap. If these excitations 
can be identified with known molecular properties, they can serve as a fingerprint for 
molecular identification. 
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Figure 2. (a) Plot of dl /dV vs Vsd demonstrating a device exhibiting conduction 
in the "weakly nonlinear" regime, (b) Plot oidl/dV vs Vsd for a device exhibiting 
a Coulomb blockade, (c) Stability diagram of a device with clear gate dependence. 
Black represent zero conductance and white represents 0.003 e 2 /h. From equation 
1, the charging energy is approximately 20 meV. The small charging energy along 
with the high number of accessible charge states indicates that the active element 
is not a single molecule of interest, but rather a metal particle, (d) Stability 
diagram of device that possesses a single Ceo molecule as the active element. 
E c > 300 meV. Dark blue represents zero conductance and red represents 0.4 
e 2 Jh. The visible excited states are consistent with vibrational excitations of the 
molecule. 



Higher-order tunneling events known as cotunneling may be observed inside the 
Coulomb blockade region of an SMT when the overall tunneling current is very small. 
Cotunneling becomes more apparent as the coupling between the dot and leads is 
enhanced, and can give rise to the non-zero current inside the blockaded region. Elastic 
cotunneling corresponds to an electron tunneling into and out of the same energy level, 
such that the molecule remains in its ground state (Fig [3j?) . Inelastic cotunneling 
occurs when an electron enters and exits the molecule through two different energy 
levels, ultimately leaving the molecule in an excited state. When eVsn = <5, where 5 is 
the energy level spacing between the ground state and the first excited state, inelastic 
cotunneling processes can occur (Figure [3f) [53,54]. Cotunneling events are only 
weakly affected by Vg and can be most clearly seen in a stability diagram mapping 
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Figure 3. (a) Energy level diagram of a single molecule transistor. A is the single 
particle level spacing, (b) Stability diagram of a model SMT displaying dl/dVsD 
(shown as relative brightness) as a function of Vgrj and Vq. A transition from 
the charge state n to the state n + 1 is shown for a molecule with n electrons, (c) 
Resonant tunneling, (d) Inelastic resonant tunneling, (c) Elastic cotunncling. (f) 
Inelastic cotunncling. (g) Cotunncling leading to the Kondo effect. 



d 2 I/dV 2 . Vibrational inelastic cotunneling is equivalent to inelastic electron tunneling 
spectroscopy [28]. 

Excited states can originate from several possible degrees of freedom, including 
excited electronic states of the system and, in molecules, internal vibrational modes. 
If the pattern of excitations is identical for multiple charge states of a given molecule 
(i.e. excitation are observed on both sides of V c at the same value of Vsd) this suggests 
that the excitations are not associated with molecular charge states and are therefore 
independent of the electronic configuration. Additionally, if excitation energies are 
equally spaced then their origin is also unlikely to be excited electronic charge states. 

As long as the exact structure of a molecule is known, its vibrational modes 
can be computed, and further calculations can demonstrate that only some of these 
modes are plausible candidates for excitation peaks observed in stability diagrams. [55] 
For instance, experimental and theoretical studies of Cqo show that the lowest energy 
internal vibrational mode energy w 33 meV and corresponds to the deformation of the 
Buckminstcrfullerene sphere into a prolate ellipsoid. [23] Similar considerations have 
been used to distinguish a known vibrational mode in C140. [29] In addition it is also 
possible to excite a center-of-mass oscillation of a molecule (commonly known as the 
bouncing ball mode) within the confinement potential that binds it to the electrode 
surface [23]. 

An additional higher-order tunneling process that is apparent in many SMT 
devices is the formation of a many-body Kondo resonance. This phenomenon has 
studied in semiconductor quantum dots [56-58] and SMTs [24-27,31,37]. Observing 
the Kondo effect requires an unpaired electron to exist in the active element. In 
the framework of the Anderson single level impurity model, the active element is 
then an effective magnetic impurity [59] with the singly occupied level's energy, e, 
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below the source/drain Fermi level (Figure QJi) ■ Because of the Coulomb interaction, 
it is classically forbidden to bring the electron out of the impurity without adding 
energy into the system. However, higher order tunneling processes can take place 
such that another electron from the source/drain electrode Fermi sea may exchange 
with the local moment. At low temperatures the coherent superposition of all possible 
cotunncling events results in the screening of the local spin; the resulting ground state 
is a correlated many-body singlet state spanning the source, impurity, and drain. This 
Kondo resonance is manifested as a conductance peak at the Fermi energy (Vsd = 0). 

The Kondo temperature, 7k, is the characteristic temperature associated with the 
formation of the Kondo singlet. The characteristic energy scale fceTk is exponentially 
dependent on V, the intrinsic width of the single particle level, by the relation [60] 

kB T K = ^l e ^(Ec-e)/TE c (3) 

outside the mixed valence regime (er < 1). The total width T = Ts + Td is dictated by 
the overlap between the single particle state and the conduction electron states of the 
source and drain. Therefore Y is in turn exponentially sensitive to the precise molecule- 
electrode configuration. This indicates both that a Kondo resonance will only be found 
in devices with relatively strong molecule-electrode coupling, and for SMTs displaying 
a Kondo resonance, the overall conductance level will be significantly increased with 
respect to similar SMTs in the Coulomb blockade regime. In a symmetric system, as 
T — ► 0, the SMT will approach its theoretical maximum conductance Go = 2e 2 /h. It 
is possible that a relatively high Kondo temperature can be associated with a device 
exhibiting a lower conductance as the total conductance is determined by the smaller 
of Tg and Td while 7k is determined by the total T. 




.40 40 10 20 30 40 50 
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Figure 4. (a) Energy level diagram of a single molecule transistor device in the 
Kondo regime, (b) Stability diagram of an SMT showing a charge state transition 
in which tunneling shifts from the Coulomb blockade regime to the Kondo regime. 
Colorbar ranges from zero conductance (dark blue) to ~ 0.5e 2 /h. (c) Traces of 
dl/dV vs Vsd displaying a Kondo resonance at several different temperatures, 
(d) Peak conductance of resonances in (c) as a function of T fit to equation 4. 
T K ~ 58 K. 

The Kondo resonance will decrease in magnitude and increase in width as 
temperature is increased. A fit to the semiempirical expression for the spin-i Kondo 
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resonance related to maximum conductance 




(4) 



where s — 0.22, or to the full width at half max 



FWHM = -^(7rfc B r) 2 + 2(fc B 7k) 2 . 



(5) 



allows extraction of the corresponding Kondo temperature, Ik- Measurements of G(T) 
or FWHM as a function of T that cannot be fit to these expressions likely indicate a 
zero-bias resonance that does not originate from a Kondo state. Zeeman splitting of 
the zero-bias resonance in an applied magnetic field can also be used to determine the 
magnetic nature of molecular state, and hence to corroborate the presence of a Kondo 
state. Finally, we always attempt to verify that the measured stability diagrams can 
be seen to transition from the Kondo regime to the Coulomb blockade regime. This 
evidence that charge state of the device can be altered with the gate together with 
the data from the tunneling thresholds of the potential SMT in the Coulomb regime 
will further attest to the nature of the device. 

The Kondo effect can be an additional "fingerprint" used to identify molecular 
states. Gate dependence of the Kondo resonance in SMTs has proven to be weaker 
than that seen in semiconductor dots [31], and "satellite" peaks associated with Kondo 
physics in the presence of molecular vibrations have been observed [26,31,61]. Also, 
since only odd occupancy charge states can exhibit the Kondo effect, the presence of a 
zero-bias resonance can be used to help identify specific charge states. Care must still 
be exercised, however: Zero-bias resonances that do not result from the Kondo effect 
may occur in some devices, while devices believed to contain metal nanoclusters [47] 
can exhibit the Kondo effect in the absence of molecules. 

4. Optical effects in electromigrated nanogaps 

Combined optical and transport experiments on electromigrated nanogaps can 
reveal a wealth of additional information beyond that available in purely electronic 
measurements. The same source/drain electrodes used to couple current to the 
nanogap are observed to act as tremendously effective plasmonic antennas, leading to 
dramatic surface enhanced Raman scattering (SERS) in the junctions [62,63]. Here we 
describe recent results in combining SERS and transport measurements, and report 
an additional effect: optical production of dc electrical currents in these molecular 
nanogap systems. These currents are consistent with optical rectification due to 
nonlinearity of the electromigrated junction's conduction, and provide a means of 
estimating the magnitude of the enhanced optical fields in the junction region. 

4-1- Surface- enhanced Raman scattering 

We have performed a series of optical measurements of Au nanogap structures 
(prepared as in Section ^ using a WITec CRM 200 Confocal Raman microscope. 
Measurements were made with a 785 nm diode laser with an incident power of 
«0.5 mW, chosen to maximize signal, minimize photodamage to assembled molecules, 
and avoid thermally driven rearrangement of the nanogap electrodes. All optical 
measurements were performed at room temperature in air. Nanogap devices were 
located in the microscope by rastering the sample beneath the microscope objective 
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to create a spatial map of the Raman response with step sizes as small as 10 nm. 
The Si substrate's strong 520 cm" 1 Raman peak can be used to map out the position 
of nanogap since the Au film of the electrodes attenuates the Si Raman response. 
Once the nanogap was located, Raman spectra were taken with 1-2 s integration time 
periods to study the dynamics of the nanogap system. 

Initial experiments examined nanogaps as a potential surface-enhanced Raman 
spectroscopy (SERS) substrate [62], with para-mcrcaptoanilinc (pMA) as the molecule 
of interest. Nanoconstrictions were placed in parallel to allow simultaneous 
electromigration of seven nanogaps at one time. Samples were characterized in the 
Raman microscope via spatial maps and time spectra of the SERS response. Prior to 
electromigration, no significant SERS response is detected anywhere on the devices. 

Following electromigration, we observe a SERS response strongly localized to 
the resulting gaps. Successive spectra measured directly over the SERS hotspot 
revealed "blinking" and spectral diffusion, phenomena often associated with singlc- 
or few-molecule Raman sensitivity. Blinking occurs when the Raman spectrum 
rapidly changes on the second time scale with the amplitudes of different modes 
changing independently of one another. Spectral shifts as large as ±20 cm" 1 were 
observed, making it difficult to directly compare SERS spectra with other published 
results. Blinking and spectral shifts are attributed to movement or rearrangement 
of the molecule relative to the metallic substrate. It is unlikely that an ensemble of 
molecules would experience the same rearrangements synchronously and thus blinking 
and wandering are expected to be observed only in situations where a few molecules 
are probed. 

More recently, we have performed simultaneous SERS and transport measure- 
ments [63], including Raman microscope observations over the center of nanogap de- 
vices during electromigration. Molecules of interest, pMA or a fluorinatcd oligimer 
(FOPE), were assembled on the Au surface prior to electromigration. A 100 x ultra- 
long working distance objective was used to allow electrical probes to be placed be- 
neath the objective to make contact with the nanogap source and drain electrodes. 
The nanogaps were migrated in situ using a computer controlled DAQ. Transport 
measurements were made by sourcing a 50-100 mV RMS sine wave at 200 Hz using 
a SRS 830 lock-in amplifcr into one electrode, with a Keithley 428 current-to-voltage 
amplifier connected at the other electrode. The ac current and its second harmonic 
were measured with lock- in amplifiers and the dc current was measured at 5.0 kHz 
using a DAQ. Simultaneously Raman spectra were acquired with a 1 s integration 
time. 

Optical measurements during electromigration provide a wealth of information 
about the plasmonic properties of nanogaps. Once the device resistance exceeds 
approximately 1 kf2, SERS can be seen. This indicates that localized plasmon modes 
responsible for the large SERS enhancements may now be excited. As the gap further 
migrates the SERS response is seen to scale logarithmically with the device resistance 
until the resistance reaches approximately 1 Mf2. In most samples the Raman response 
and conduction of the nanogap become decoupled at this point with the conduction 
typically changing little while uncorrelated Raman blinking occurs. A more extensive 
discussion of the connection between plasmonic modes and interelectrode conductance 
is presented elsewhere [63]. 

In about 11% of 190 devices, however, the Raman response and conduction show 
very strong temporal correlations. The 11% yield is is quantitatively consistent with 
yields of gateable SMTs as mentioned in Section [3] A typical correlated SERS time 
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Figure 5. Waterfall plot of Raman spectrum (1 s integration) and correlated 
conductance measurement for a single FOPE molecule shown at bottom. The 
Raman mode observed between 1950 cm" 1 to 2122 cm" 1 is believed to be for the 
same 2122 cm" 1 mode associated with the C=C stretch of the FOPE molecule. 
The large spectral shifts observed for this mode are attributed to interactions 
between the molecule and its nanogap environment. Clear correlations between 
the Raman structure and conductance can be seen. In particular during region B 
and for part of region E the Raman spectrum is observed to disappear while the 
conductance drops to zero. Between different regions distinct changes in spectrum 
are observed during clear changes in conductivity. 

spectrum and conductance measurement for a FOPE device are presented in figure 
[5] The temporal correlations between SERS and conduction are clear. In region A 
we observe a stable Raman spectrum and small conductance changes which appear to 
be correlated with changes in the Raman mode at 1980 cm -1 . After a short spike in 
conductance the conductance and Raman disappear at B. In section C the conductance 
and Raman spectrum return. This time the conductance is 6x larger than in A and 
the Raman modes are different. In particular the mode previously seen at 1980 cm" 1 
is now at 1933 cm -1 . In region D we see the conductance drop to levels closer to those 
seen at A and the 1933 cm -1 mode from C shifts to 2038 cm" 1 . During D we see that 
the conductance momentarily returns to the value seen in section C correlated with a 
shift in the 2038 cm" 1 mode back to 1933 cm -1 . In section E we see another shift in 
position of the 1933 cm" 1 mode to 2098 cm" 1 which slowly shifts up to 2122 cm" 1 . 
During E the Raman spectrum is seem to once again disappear correlated with drop in 
conductivity briefly and the return simultaneously. Finally in F we see the spectrum 
disappear a final time again correlated with a drop in conductivity. 

In the bulk Raman spectrum the FOPE molecule only shows one mode above 
1700 cm -1 , the C=C stretch mode at 2228 cm -1 associated with the triple bond 
connecting the two phenylene rings. It is likely that the mode at 2122 cm -1 is a 
manifestation to the 2228 cm" 1 mode. This mode frequently appears when studying 
the FOPE devices and is absent in control experiments with pMA or contaminants 
from the air. The large shift in wave numbers (over 100 cm -1 between the bulk and 
normal and almost 300 cm -1 for the greatest shift) indicates significant interactions 
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between the molecule and either the substrate or other adsorbates, and is cause for 
further investigation. 

Recalling that conduction in nanogaps is dominated by approximately a single 
molecular volume, the observed correlations between conductance and Raman 
measurements strongly indicate that the nanogaps have single molecule Raman 
sensitivity. It is then possible to confirm that electronic transport is taking place 
through the molecule of interest, via the characteristic Raman spectrum. Data sets 
such as Figure [5] contain implicitly an enormous amount of information about the 
configuration of the molecule in the junction, and should be amenable to comparisons 
with theoretical calculations of the optical properties of the molecule/electrode region. 
These combined optical and transport measurements open many possible paths of 
exploration. 

4-2. Optically Induced Transport 

In addition to SERS, we also observe significant dc currents in electromigrated 
nanogaps under illumination. The mechanism of this optically induced transport 
may be probed by studying the dependence of the dc current on the incident 
optical intensity and the measured low frequency transport properties of the junction. 
Resonant optical effects [64,66] and photon-assisted tunneling [67] are potential sources 
of dc optically-driven currents in molecular junctions. However, we find that the 
optically-driven dc currents are relatively independent of molecule type, with similar 
data sets collected in devices using, for example, pMA molecules, FOPE molecules, 
and adsorbed atmospheric contaminants. This strongly suggests that the mechanism 
behind these optically-driven currents is a general feature of the electromigrated 
junction structure, rather than tied to specific molecular features. 

One mechanism that is consistent with our observations is optical rectification 
due to the nonlinearity of the source/drain tunneling characteristics. This effect has 
long been considered in STM experiments [68], though its unambiguous observation 
has been very challenging [69]. The rectified current originates from the interaction 
of an ac excitation with a nonlinear circuit element. In the limit of a small bias V 
the current can be approximated via a Taylor series. In particular for an oscillating 
potential V = Vo + V ac cos(ojt): 

^ ac 2 C0S 2 (^) + ...(6) 

Vo 



V; c 2 (cos(2wt)+l) + . . . 

Vo 

The conduction nonlinearity leads to a second-harmonic ac signal as well as an 
additional dc current, both linearly proportional to d 2 I/dV 2 , which will depend on 
the device geometry and conduction through the molecule. Additionally the optically 
rectified current will depend linearly on the incident laser intensity. 

Note that optical rectification in nanogap devices would allow an experimental 
estimate of the enhanced optical field. One can measure d 2 I/dV 2 using a low 
frequency (e.g., 200 Hz) ac signal. If this nonlinearity results from tunneling and 
the tunneling timescale is fast compared to an optical cycle, then one can use the 
dc optically-driven dc current to infer V op t, the optical- frequency potential difference 



W=W) + (§)^cos W+ i(|4) 

Applying a trigonometric identity, 

, x ( dl \ , s 1 / d 2 I \ 
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Figure 6. Schematic of optical/electronic measurement. A 100 mV RMS ac 
signal is applied by a lock-in into the source electrode at frequency lu. The ac 
current is measured at the drain by a current-to- voltage amplifier. The dc current 
is measured from the amplifier by a DAQ at 5 kHz. The ac current is measured 
at u) and 2ui using lock-in amplifiers. Raman spectra are acquired using 1 s 
integrations at an incident power of 0.5 mW at 785 nm. For later experiments 
the equipment enclosed by the dashed lines was added to allow the laser to be 
chopped at frequency ujo- <^ and ljo have been chosen such that they and 2u> are 
at least 50 Hz apart. The optically rectified current is measured via a lock-in at 
the chopping frequency. 



across the source/drain electrodes at the point of tunneling. With an estimate of the 
source/drain gap distance from dl/dV, an estimate of the local plasmon-enhanccd 
electric field is then possible. 

Using the measurement scheme shown in Figure|6]wc measure the first and second 
harmonics of the low-frequency current using a lock-in amplifier, as well as the dc 
current. In the absence of any optical effect. ([7]) implies that the rectified current 
due to the low-frequency drive will be exactly equal to the lock-in current measured 
at the second harmonic of the source frequency, as we observe. Under illumination, 
plots of the dc current vs. the measured low-frequency ac current at 2uj for fixed 
input amplitude of 100 mV show that the dc current can exceed the low frequency dc 
contribution by nearly a factor of two. This implies a second source of dc current that 
scales linearly with d 2 I/dV 2 and vanishes when the illumination is blocked. 

To determine the particular optical mechanism at work, we measured the 
dependence of the dc current on optical power. We chop the incoming laser and 
measure the current component at the chopping frequency as shown in the Figure 
[6] Figure [7] shows a power curve measured at a few intensities, showing that the 
optically induced dc current depends linearly on the incident intensity, consistent 
with the optical rectification mechanism. 

Assuming that optical rectification is at work we can estimate the enhancement 
of the electric field in the nanogap, as described above. By comparing the measured 
dc current to the low frequency ac current component at 2lo at constant laser power 
we can determine how much of the rectified current is a result of optical rectification. 
Figure [7)3 shows a representative curve. The dc current is well fit by a line of slope - 
1.91. A slope of ±1.0 would have indicated that all the dc current is due to rectification 
of the applied low frequency signal. The sign of the slope varies from device to device 
depending on the sign of d 2 I/dV 2 and thus the direction of the rectified current. It 
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Figure 7. a) Rectified current as a function of optical intensity. The measurement 
was performed on a nanogap assembled with pMA and migrated to a resistance 
of ~10 MSI. Error bars are one standard deviation from the mean value. The 
data is well approximated with a linear fit consistent with the proposed model 
for optical rectification, b) dc current measured as a function of the ac current 
component at 2lo (source frequency). Each point averaged over 2 seconds with 
the error bars indicating one standard deviation in the dc measurement over that 
time. The data is reasonably represented by a linear fit with slope —1.91 ± 0.14 
with a r 2 value of 0.883. This slope magnitude in the presence of illumination 
implies significant optical rectification. 

should be noted that both the low frequency and optical frequency parts should be 
rectified in the same direction. A slope of -1.91 indicates that (V^pticai/Vo) 2 = 0.91. 
For this particular measurement Vb=100 mV yielding V op ti C ai=95 mV. Assuming a gap 
separation of 1 nm we can get the approximate optical field strength across the gap 
is ?s 1 x 10 8 V/m. The incident unenhanccd optical field is w 2 x 10 5 V/m, yielding a 
field enhancement of approximately 500, and a total Raman enhancement of roughly 
6 x 10 10 , consistent with the predications for the necessary enhancement to observe 
single molecule SERS. 

Unfortunately at room temperature and under higher intensity laser powers 
the nanogap stability is poor, making it difficult to perform a more in-depth 
analysis. Currently we do not have the capability to perform low temperature optical 
measurements; however, it is possible to examine these same rectification mechanisms 
at radio frequencies in a low temperature probe station. Optical rectification has been 
successfully demonstrated at microwave frequencies in STM measurements [69]. A 
measurement scheme analogous to that present in Figure [6] was used, with d 2 I/dV 2 
found by a low frequency measurement and the RF contribution to the rectified current 
measured by chopping the incident RF. 

Devices were prepared by assembling pMA on the Au surface prior to the 
experiment. Measurements were made at 80 K in a vacuum probe station after 
electromigration. A 20 mV RMS low frequency sine wave was sourced by one lock-in 
amplifier and fed into a bias tee along with a 0-20 mV RMS RF signal. The RF signal 
was amplitude modulated from to 100% a frequency uj (typically 157 Hz). The 
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Figure 8. (a) RF Rectification vs \d 2 I/8V 2 V n¥ measured at 10 MHz. The 
rectified current scales linearly as expected by our model. After factoring in 
the impedance mismatch between the RF function generator and the nanogap 
structure the data is fit well by a line of slope=l. This experimentally confirms 
equation [7] as the correct mechanism for rectification in nanogaps. (b) d^I/dV 2 
as a function of the applied RF amplitude. As the RF amplitude is increased the 
average nonlincarity increases in a linear fashion. 

combined signal was sent into one electrode of the nanogap. The current flow through 
the device was measured at the other electrode using a current amplifier with 100 kHz 
low pass filter enabled to remove any RF component from the measured current. 
The output of the current amplifier was measured using three lock-in amplifiers 
measuring at uj, 2cj, and ujq (lock-in source frequency, second harmonic, and RF 
chopping frequency) and the dc current was measured directly at the amplifier. A 
representative power curve plotting j(d 2 I/dV 2 ) Vrp versus I UQ is plotted in figured] 
for an RF frequency of 10 MHz. A clear linear dependence is observed. The sloped of 
4.48 is a result of an impedance mismatch between the RF function generator and the 
nanogap device. After accounting for this mismatch the slope is approximately 1 as 
expected for the rectification mechanism. From this series of measurements we infer 
that d 2 I/dV 2 is frequency independent up to 1 GHz. A priori this does not imply 
that d 2 I/dV 2 is frequency independent up to optical frequencies, though the results 
seen in figure [7] are consistent with this rectification mechanism beyond ~ 10 14 Hz. 

5. Conclusions and prospects 

Electromigrated gaps have proven to be an enormously useful tool in probing electronic 
properties in single molecules, though the lack of imaging techniques has meant 
that great care must be exercised in interpreting transport data. Pure electronic 
transport can contain signatures that are distinct to the molecules being probed 
(charging energies, vibrational resonances, Kondo physics). The recent observation 
that electromigrated junctions are highly effective optical antennas has great potential. 
Simultaneous measurements of transport and Raman spectra in single molecules are 
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now possible, allowing the vibrational fingerprinting of the active electronic element 
in the junction. Detailed information about molecular orientation and bonding within 
the junction may be inferrable from the Raman data. Optical measurements also 
bring additional electronic transport mechanisms into play, as optical rectification 
measurements demonstrate. The promise of combined single-molecule electronic and 
optical characterization suggests that electromigrated molecular junctions have a 
bright future as tools for physics and physical chemistry at the molecular scale. 
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